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Using standard definitions of chaos (as positive Kolmogorov- Sinai entropy) and diffusion (that 
multiple time distribution functions are Gaussian), we show numerically that both chaotic and 
^— ^ ", nonchaotic systems exhibit diffusion, and hence that there is no direct logical connection between 

the two properties. This extends a previous result for two time distribution functions. 
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I. INTRODUCTION 



It is generally surmised that microscopic chaos is a necessary condition for diffusive behavior of a system. Microscopic 
chaos here means a positive Kolmogorov-Sinai (KS) entropy, or equivalently, at least one positive Lyapunov exponent. 
While this seemed to be confirmed by an ingenious experiment with a Brownian particle [Q, we showed that the 
same behavior on which the microscopic chaoticity of the Brownian system was established also held for a manifestly 
nonchaotic system with zero KS entropy so that both systems exhibited diffusion. [] In these works we used the 
f*^ \ mean square displacement (Eq. |^ below) , a two time distribution function, as our definition of diffusion. 

An anonymous referee of a Ref. || , alerted us to a definition of diffusion (Eq. |l|) which involved multitime distribution 
functions and the question arose whether those could distinguish between chaotic and nonchaotic models, in that all 
our chaotic models and none of our nonchaotic models would satisfy this definition. In this note we show that this 
behavior does not occur: no difference has been found between the diffusive behavior of microscopically chaotic and 
ii . nonchaotic models even in the multitime distribution functions we studied. 

Our models, a subset of those of Ref. ||], consist of a point particle moving with unit speed in two dimensions, un- 
dergoing specular reflections (collisions) with fixed non-overlapping scatterers. The shapes, positions, and orientations 
of the scatterers in the various models are defined in the next section, and can be classified as Lorentz or (modified) 
Ehrenfest models. For each fixed configuration of scatterers, the particle's position x is determined numerically as a 
function of the time t for random initial positions x(0) (which are uniformly distributed over the plane, except that 
' they cannot lie inside a scatterer) and random initial velocity directions (also uniformly distributed) of the particle. 
i-H , The probabilities below are defined with respect to this distribution of initial conditions. 

The definition of diffusion proposed by the referee and used in this paper is stated in terms of the probability of 
the particle initially at the origin visiting n regions at n different specified times. Explicitly, a system is said to be 
diffusive, if, for all n — 1, 2 . . ., the n + 1 time (including t = 0) probability distributions V satisfy the relation j7j: 

• i— ( ' 
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where < t\ < i 2 < . . . < t n are arbitrary times and the T>k (k = 1, 2 . . . n) are arbitrary subsets of the plane. p(x, t) 
is the Gaussian Green's function of the (generally anisotropic) diffusion equation in two dimensions, 

e -ti 2 /4B„„t-u 2 /4C„i 

^ ' 47rV D uu D vv t 

and determines the probability density of a displacement (Au, Xv) in time \ 2 t in the limit A — ► oo. The scaling of the 
position and time by A is based on the fact that for Gaussian diffusion the mean square displacement of the particle 



1 There are at least two integrable models (hence, with zero Lyapunov exponents) that have been shown to exhibit diffusion 
as defined by Eq. [j], a one dimensional hard rod model and a heavy particle in a lattice of light particles || . In both cases 
there are an infinite number of moving particles, so we expect that the KS entropy is infinite, as in the infinite ideal gas Q. 
Here we have only one moving particle, and we argued in Ref. |3[ that the KS entropy of our nonchaotic models is zero. 
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is proportional to the time, and is a macroscopic limit. Here, u and v are distances measured in the eigendirections 
of the diffusion tensor, while D uu and D vv are the corresponding eigenvalues. In the isotropic case D uu — D vv = D 
and the choice of orthogonal directions (it, v) is arbitrary. 

The components of the diffusion tensor in an arbitrary orthonormal basis {e^} (so that x = y\ Xjej) can be related 
to the mean square displacement — 1,2): 



Dij = i im _ ^±21 (3 ) 



t^oo 2t 

where Ax = x(i) — x(0). This expression involves the usual 2-time distribution function, that is, n = 1 in Eq. (p]). 
The mean square displacement is discussed in Ref . || , where we found numerically that both chaotic and nonchaotic 
models exhibited Gaussian diffusion. It is conceivable that a system may have Gaussian 2-time distribution functions 
but non-Gaussian multi-time distribution functions, i.e., for larger values of n. The results of this note show that also 
for larger n, both chaotic and nonchaotic models appear diffusive within the accuracy of the numerical methods. 



II. MODELS 

The models we consider are the same as some of those discussed in detail in Ref. As mentioned above, they all 
contain a point particle colliding with fixed non-overlapping scatterers in two dimensions. The density is such that 
the scatterers cover exactly half of the available area. 

The microscopically chaotic models are Lorentz models, consisting of circular scatterers distributed in position 
either periodically or randomly over the plane, such that they do not overlap. They are microscopically chaotic as 
they have a positive KS entropy. See Fig. |] 

The microscopically nonchaotic models are Ehrcnfcst-like models, consisting of square scatterers distributed in 
position either periodically or randomly over the plane, such that they do not overlap. They are microscopically 
nonchaotic, with zero KS entropy. Unlike the original Ehrenfest model, the squares in this paper are always oriented 
at random. See Fig. [j]. 
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FIG. 1. The periodic Lorentz model (left) and periodic modified Ehrenfest model (right). 

In both the Lorentz and Ehrenfest periodic models, the elementary cell is square and contains four scatterers. The 
random models have no preferred directions, so the diffusion is isotropic. In contrast, the periodic models have a 
directionality defined by the elementary cell, and the diffusion is anisotropic as in Eq. (g). 
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2 Specifically: the periodic circles were denoted LP4 in ||, the randomly positioned circles Loo, the periodic squares RP4 and 
the randomly positioned squares Roo. Here, the L stands for the Lorentz gas, R for the random orientations of the squares, P 
for periodic, 4 for the linear dimension of the unit cell (there are also 4 scatterers per unit cell), and oo for the infinite "unit 
cell". 
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We investigated Eq. ([j]) for equally spaced times and two choices of the subsets T>k of the plane: 
A The T>k (k < n) are all balls of radius e = l/\/T0 centered on the origin for n = 1 . . .7. 



B n = 2, where and are two balls of radius e = 1/VTO centered on randomly chosen points xj?; i = 1 ... 7 
identifies the different members of our sample of seven configurations of the T> k , see Table § 
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TABLE I. Positions defining the centers of the Tj£ for case B, and corresponding symbols in Fig. ^ below. 



In both cases the scale factor varied as A = 1 . . . 100, expanding the times by A 2 and the subsets T) k by A as 
required by Eq. (Q). In case A this means that as A increases the T> k become larger circles centered on the origin, 
while in case B the T> k become larger but also move off to infinity. The time steps were a single time unit, so that 
tfc = k (k — 1 . . . n), and the maximum length of a trajectory was X 2 t n = 70000 for A = 100 and n = 7. We took 
1000 trajectories of length 10 6 with different initial conditions. In order to get reasonable statistics without excessive 
computer time we considered of order 10 6 segments of length 7 x 10 4 along each trajectory, shifting the initial state 
on the trajectory by one time unit. The sample size was thus of order 10 9 . Because these trajectory segments overlap 
and hence are not independent, the statistical error is difficult to estimate a priori; we find below that it is greatest for 
the periodic square model, which is not surprising given that it is expected to have the largest amount of correlation. 

In A we are testing for the probability for the particle to return n times at regular intervals to near its starting 
point. The motivation is that while in nonchaotic systems, a trajectory can follow an almost periodic trajectory for 
a long time, leading to a power law decay of recurrences, in chaotic systems the decay of recurrences is exponential. 
The power law for nonchaotic systems could lead to correlations over long times, inconsistent with the uncorrelated 
diffusion of Eq. (Fy). The integrand in Eq.fiLj) is a product, so the successive time steps are conditionally independent. 



In B we consider alternative choices for the T>k in the (we think) unlikely event that the rather special case A has 
different long time behavior than the general case. Case B does not extract almost periodic orbits but it may reveal 
other ways in which correlations might be manifest at long times. 

III. RESULTS 

We measured the following quantities for the cases A and B discussed above: 

1. The probability V{n\ A; {V k }) for various A, n and V k as discussed above. We often omit some of the arguments 
for simplicity. 



3 Case A is similar to the "almost periodic recurrence" method of Ref. M, which showed a difference between chaotic and 
nonchaotic models by looking at the probability for return to near the starting point many times at regular intervals. The 
difference is that in Ref. Q the size of the regions remained fixed as the time of return increased, while in this paper the T> k 
and the times are scaled by factors A and A 2 , respectively. 

4 It would be somewhat impractical to go beyond n — 2 in case B: (a) the integral on the RHS would be of increasingly high 
dimension and hence hard to evaluate, and (b) the probabilities would be much lower so that the statistics would be very poor, 
see Fig. below. 
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2. For case A, the conditional probability V c (n + l\n; A; = Vin + 1)/V(ri) for the particle being in XD n+ i 

given that it was in all the previous XDk for fc = 1 . . . n. Again, we often omit arguments for simplicity. 

Eq. (|j) actually makes two testable assertions: Does the limit on the left hand side (LHS) exist, and if so, is it 
equal to the right hand side (RHS)? The above measured quantities give information about the LHS and we are 
first interested in the large A limit of these quantities; for the RHS we obtain the diffusion tensor from Eq. (|J) (see 
Table |n]) and perform the integration numerically. This has only been carried out up to n = 2 (which is a four-fold 
integration). However, for larger n and case A, the right hand side is (to a very good approximation) multiplied 
by a constant factor for each additional repetition, ^| so the prediction is that the conditional probabilities become 
independent of n for large A and n. In fact our result show this behavior already at the smallest values of n. 
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TABLE II. Components of the diffusion tensor computed in original coordinates (x,y) from Eq. (|3J), and 
then diagonalized to give the (u, v) values. Note that the randomly placed and oriented squares as well as 
the randomly placed circles are isotropic, with D xy = and D xx = D yy within errors of about 0.001, while 
the periodic models are anisotropic and satisfy neither of these conditions. 



The results are shown in Figs. UH and Table [II. For all n and {T>k} the V{\) approach a constant (apart from 
statistical error) at large A, indicating convergence of the limit in Eq. ([|), see Fig. g for case A and Fig. ^ for case B. 
For case A, all the conditional probabilities approach the same constant indicating conditional independence of the 
multi-time probabilities as predicted by diffusion, see Fig. |^; an equivalent statement is that on a logarithmic scale, 
all the curves in Fig. || become equidistant at large A. We have verified numerically that in both cases A and B the 
values of the constants as A — * oo are consistent with their values as given on the RHS of Eq. (0), thus confir mine the 
diffusive nature of all these models. 
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TABLE III. Conditional probabilities Vc(n + l\n). The second, "theoretical" column gives the values 
calculated from numerical integration of the RHS of Eq. (|l]) . The other columns, derived from Fig. ^| by 
averaging over A from 50 to 100, agree with the values from column 2, and do not depend on n, consistent 
with conditional independence in the underlying probability distribution. The differences in the periodic 
square model are almost certainly due to statistical errors, as in Figs. B, 0. 



5 The approximation can perhaps be understood as follows. If we perform all but the last x n integration, we are left with an 
integral over a function, say f n (x. n ). At the n + 1 stage, we have f n +i as a linear operator C (integral over a Gaussian kernel) 
acting on f n . Assuming £ is well behaved, the /„, approach a limiting form /Aq as n — > oo where < Ao < 1 is the leading 
eigenvalue of C and / is the leading eigenfunction. Performing the last integration, we find that successive integrals differ by a 
factor of Ao, leading to a vanishing probability for n — *• oo. 
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FIG. 2. Probabilities V(X) as a function of A for case A, with different curves giving n — 1 . . .7. Note 
that the probability becomes independent of A indicating the existence of the limit on the LHS of Eq. (|l|). 
The equal spacing between subsequent values of n indicates a fixed ratio of successive recurrences since the 
scale is logarithmic. This is consistent with the uncorrelated diffusion as given by the RHS of Eq. (|l|). The 
statistical error is greatest when the probability is low, and when the correlations in the data are greatest, 
that is, for the periodic square model. 
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FIG. 3. Conditional probabilities Vein + l|n; A) as a function of A for case A, and n = 1 ... 4, derived 
from Fig. ^. The equal spacing in Fig. |^ now causes the graphs corresponding to different values of n to 
coincide. Again, statistical error is greatest with low probability (larger n) and with greater correlations in 
the data (the periodic square model). 
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FIG. 4. Probabilities V{X) as a function of A for case B. The different symbols correspond to different 
2?jj* ' (i — f . . . 7) as given in Table [E| The dashed or dotted horizontal lines give the predicted values from 
the RHS of Eq. (Q), and are not visible when they run through the data. The order of the curves differ due 
to the anisotropy of the diffusion in the periodic models. 



IV. DISCUSSION 



1. The above mentioned results suggest that there is no difference in the diffusive behavior of chaotic (Lorentz) 
and non-chaotic (Ehrenfest) models. Both exhibit Gaussian diffusion and their only difference is the value of 
the diffusion coefficients. 

2. The case A, in particular the larger probabilities for recurrence at small A as seen in Fig. especially suggests 
the presence of "small" periodic orbits, which favor return to previously visited regions for the nonchaotic wind- 



7 



tree models since they exhibit only algebraic 1/t decay (see Ref. ||). For the chaotic Lorentz models, periodic 
orbits disappear exponentially fast in t and therefore more returns are less likely. Nevertheless this difference in 
periodic orbit stability, which indicates a more subtle chaotic difference than the corresponding KS entropies in 
these models, do not appear to lead to a qualitatively different, non-Gaussian, diffusion process. 

3. The influence of small almost periodic orbits is clearly local and incoherent from region to region, ft is therefore 
constantly wiped out and does not influence the global diffusion process at all, except for example in the values 
of the Dij, which should decrease by a strong influence of almost periodic orbits. 

4. Our results here appear to differ from the behavior observed in open (that is, finite with absorbing boundaries) 
Lorentz and Ehrenfest models in Ref. ||. There, we studied the number of particles with different initial 
conditions remaining in the system as a function of time. Only the Lorentz model exhibited the exponential 
decay of particles characteristic of the diffusion equation; the Ehrenfest model corresponding to our squares 
exhibited l/t decay due to trapping of particles around periodic orbits. This is not in contradiction to our 
results here (for which the nonchaotic model also exhibits diffusion) since the limits are different: With absorbing 
boundaries, the limit t — > oo occurs with the linear dimension of the observed region (the whole system) L fixed; 
in this note the linear dimension of the observed regions XD is scaled as t 1 ^ 2 , which is more natural for diffusion. 
Thus, the identification of a model with diffusion depends not only on the microscopic dynamics, but also on 
the manner in which the limit is taken, ft would be interesting to investigate these relationships further. 

5. To conclude, our results indicate that at least for the models we consider, subtle correlations, whether associated 
with almost periodic orbits in small clusters of scatterers or otherwise, are insufficient to modify the long time 
Gaussian form of the diffusion process when the natural macroscopic limit is taken (such that lengths are scaled 
with the squareroot of times). In so far as these correlations are determined by the microscopic chaotic nature 
of the system, the chaoticity does not seem to show up in the macroscopic diffusive behavior. This behavior 
is consistent with the results of the Brownian motion experiment discussed in Ref. [Q, and the argument || 
that the long measuring times of l/60s means that the results could come equally from chaotic or nonchaotic 
microscopic dynamics, and that any experimental determination of microscopic chaoticity would necessarily 
involve measurements at microscopic time scales. 
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